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We have constructed a theoretical framework of the biexciton-resonant hyperparametric scattering 
for the pursuit of high-power and high-quality generation of entangled photon pairs. Our frame- 
work is applicable to nano-to-bulk crossover regime where the center-of-mass motion of excitons and 
biexcitons is confined. Material surroundings and the polarization correlation of generated photons 
can be considered. We have analyzed the entangled-photon generation from CuCl film, by which 
ultraviolet entangled-photon pairs are generated, and from dielectric microcavity embedding a CuCl 
layer. We have revealed that in the nano-to-bulk crossover regime we generally get a high perfor- 
mance from the viewpoint of statistical accuracy, and the generation efficiency can be enhanced by 
the optical cavity with maintaining the high performance. The nano-to-bulk crossover regime has 
a variety of degrees of freedom to tune the entangled-photon generation, and the scattering spectra 
explicitly reflect quantized exciton-photon coupled modes in the finite structure. 

PACS numbers: 42.65.Lm, 42.50.Nn, 71.35.-y, 71.36. +c 



I. INTRODUCTION 

Entangled photon pairs have been discussed in relation 
with the Einstein-Padolsky-Rosen (EPR) paradox^ and 
nowadays they play an important role in quantum infor- 
mation technologies. The pursuit of their high-quality 
and high-efficiency generation is a fascinating subject 
in the fields of quantum optics and solid-state physics. 
In addition to the standard generation method by para- 
metric down-conversion (PDC) in second-order nonlinear 
crystals^ the generation scheme using a semiconductor 
quantum dot£~— attracts much attention, because purely 
a single pair of entangled photons is created in princi- 
ple, and it can be a deterministic source of entangled 
pairs. Recently, the generation efficiency is highly en- 
hanced by implementing an optical cavity structure with 
distributed Bragg reflectors (DBRs)£ and by a molecule 
of micropillarsjS. Further, the emission by electric injec- 
tion has been reportedi 7 - On the other hand, the develop- 
ment of entangled photons as an excitation light source 
is of growing importance for the next-generation tech- 
nologies of fabrication and chemical reaction^ For this 
purpose, high-power and high-quality entangled-photon 
beams are absolutely necessary, and this high-power but 
probabilistic generation is another direction of research in 
addition to the deterministic generation by single quan- 
tum dots. 

In the process of PDC^^ an incident photon with 
frequency Wj n and wavenumber k m splits into two pho- 
tons (wi,fci) and (w2,fc2) with satisfying the conserva- 
tion of energy u)- m = lj\ + L02 and of wavevector fcj n = 
k\ + &2- This second-order nonlinear process creates 
polarization-correlated entangled-photon pairs in nonlin- 
ear optical crystals with a birefringence. On the other 
hand, Savasta et al^ suggested and Edamatsu et al.— 
experimentally demonstrated that ultraviolet entangled- 
photon pairs are generated by the biexciton-resonant hy- 
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FIG. 1: The biexciton-resonant hyperparametric scattering 
(RHPS) is depicted on dispersion curves of biexciton, exciton- 
polariton, and longitudinal exciton. Biexcitons are reso- 
nantly created by two-photon absorption, and entangled- 
photon pairs are emitted when the biexcitons decay into the 
lower exciton-polariton branch. This emission appears in 
scattering spectrum as two peaks called LEP and HEP (lower 
and higher energy polaritons) as seen in Fig. [3] Due to the 
conservation of energy and wavevector, the positions of the 
two peaks depend on scattering angle , 10 ' 11 and the entangled 
two photons are emitted symmetrical about the pump beam 
as shown in the inset. Two additional peaks called Mt and 
Ml in scattering spectra originate from the biexciton decay 
into transverse and longitudinal exciton levels, respectively. 



pcrparametric scattering (RHPS) in CuCl (see Fig. [TJ. 
The RHPS is a third-order nonlinear process, in which 
two incident photons resonantly creates a biexciton (ex- 
citonic molecule) with (2wj n , 2fcj n ) and it spontaneously 
collapses into a photon pair satisfying 2uj- ln = lo\ +W2 and 
Since the lowest level of biexcitons in 
CuCl, which was resonantly excited in the experiment, 
has zero angular momentum^ the emitted pair consists 
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of left- and right-circularly polarized photons conserving 
the total angular momentum. Owing to the two possi- 
ble decay paths involving exciton-polaritons, the emitted 
photons are polarization-entangled. 

The generation efficiency of RHPS is much higher than 
that of PDC, because of the giant oscillator strength 
of the two-photon absorption involving the biexcitonAi 
However, in the first experiment^ a part of observed 
pairs has no entanglement, and this noise was subtracted 
in the estimation of entanglement of the generated pairs. 
As indicated by Oohata et al.ji 3 - the main contribution 
of the unentangled pairs is an accidental collapse of two 
biexcitons, and this problem has been successfully sup- 
pressed by using high-repetition and weak-power laser 
pulses, because the number of unentangled pairs (noise) 
is increased by Jj n 4 for increasing the pumping power I ln 
while the number of entangled pairs (signal) is propor- 
tional to Ii n 2 - However, this fundamental trade-off prob- 
lem between signal intensity and S/N ratio should be 
resolved from the improvement of material structures 1 ^ 
in addition to the improvement of pumping condition of 
Ref. [T3I While one solution is using a single quantum 
dot as a deterministic source,— ~— for the pursuit of high- 
power generation there is a proposal of using an opti- 
cal cavity embedding an excitonic quantum well for the 
improvement of generation efficiency 15 ' 16 Furthermore, 
owing to the rapid radiative decay by the exciton super- 
radiance (enhancement of interaction volume between ex- 
citons and photons) ) 17 ' 18 we have theoretically revealed 
that the trade-off problem can be resolved simultaneously 
realizing a high generation efficiency by using an opti- 
cal cavity embedding an excitonic layer in nano-to-bulk 
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crossover regime i« 

In a microcrystal, such as quantum dot and quantum 
well, smaller than the Bohr radius of excitons, the elec- 
tron and hole are individually confined in the crystal, 
and the relative motion of excitons and also the binding 
energy are strongly modified from those in bulk crystal. 
When the crystal size is larger than the exciton Bohr ra- 
dius but small enough compared to the light wavelength, 
the center-of-mass motion of excitons are confined, and 
the center-of-mass kinetic energy is quantized ! 19 ' 20 When 
the crystal size is comparable or a few times larger than 
the wavelength (nano-to-bulk crossover regime) , the sys- 
tem is characterized by exciton-photon coupled modes 
with peculiar resonance energy and radiative life time, 
and the coupled modes are gradually reduced to bulk 
polaritons with increasing the crystal sizei 18 ' 21 " — In this 
crossover regime, the system shows a variety of optical re- 
sponses compared to bulk materials and also to quantum 
dots due to the center-of-mass confinement of excitons 
and the spatially resonant coupling with electromagnetic 
fields. Actually, owing to the recent development of nano- 
scale fabrication, anomalous nonlinear optical processes 
have been reported in semiconductor nano-structures and 
in the nano-to-bulk crossover rcgimei 17 ' 26 " — Further, it 
shows a rapid radiative decay rate of excitons on the 
order of 100 fs due to the exciton superradianceJi Con- 



cerning the entangled photon generation, while the per- 
formance of PDC method is almost governed by the 
choice of nonlinear materials and its thickness, the RHPS 
method significantly depends on the quantum states of 
excitons and biexcitons, because it is a resonant process 
involving the elementary excitations. In the nano-to- 
bulk crossover regime, the generation of entangled pho- 
ton pairs by RHPS can be significantly modified with 
respect to frequencies, angles, polarizations, and phase 
difference of the generated entangled state as discussed 
in our previous letter In the present long paper, we will 
show the detailed theoretical framework for the investiga- 
tion of the entangled-photon generation in nano-to-bulk 
crossover regime with multilayer structures, especially an 
excitonic layer embedded in DBRs. 

We explain our theoretical framework in Sec. HH and 
show in detail the way to calculate the one-photon scat- 
tering intensity and the two-photon coincidence intensity 
of RHPS in the case of multilayer structure in Sec. IIIII 
The calculation results are shown in Sec. IIV1 and the 
discussion is summarized in Sec. |V] 



II. THEORETICAL FRAMEWORK 

The emission spectra from Bosc-Einstcin condensation 
of biexcitons were calculated by Inoue and Hanamuraj^ 7 - 
and they also showed the relation between energies and 
scattering angles of two peaks called LEP and HEP 
(lower and higher energy polaritons). Later, Hanamura 
and Takagahara 3 -^ calculated line shapes of the so-called 
Mt and Ml peaks, which are emitted by the relaxations 
of biexcitons to transverse and longitudinal excitons, re- 
spectively. The entanglement of the scattered photons by 
RHPS was first pointed out by Savasta et al.ji 2 - and their 
theoretical framework^ 9 - is based on the quantum elec- 
trodynamics (QED) theory for dispersive and absorbing 
medi a 40 ' 41 and on the exciton-exciton correlation func- 
tions calculated from first principles ! 42 ' 43 

In the present paper, in order to correctly treat the 
center-of-mass confinement of excitons, we start from the 
QED theory of excitons^ which simultaneously solves 
the equation of motion of excitons and of electromag- 
netic fields with inheriting the concepts of the above QED 
theorie a 40 ' 41 and of the semiclassical nonlocal theor y 30 ' 45 
(or the so-called ABC-free theory^-). It is well known 
that the center-of-mass motion of excitons raises more 
than one propagating modes of exciton-polaritons in their 
band gap frequency, and the RHPS process has been 
used to observe the dispersion of polariton o 11 ' 37 ' 47 " — and 
also to measure the translational masses of excitons and 
biexcitons i 1 1 1 50 ^ Moreover, optical responses explicitly 
reflects the confinement of center-of-mass motion of ex- 
citons in nano-structured materials and also in the nano- 
to-bulk crossover regime jlZ r 20 1 26 ~ .36 wn i cn we investigate 
in the present paper. 

Concerning the treatment of biexcitons, we suppose 
the excitons as pure bosons and consider an exciton- 
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exciton interaction leading to the creation of biexcitons. 
However, instead of the detailed treatment in the the- 
ory of Savasta et al.,— we simply assume the relative 
motion of the lowest level of biexcitons with some pa- 
rameters measured in experiments 1 1 1 1 53 1 54 and the coeffi- 
cients of the exciton-exciton interaction is replaced by the 
wavefunction and the binding energy of biexcitons. This 
treatment is very simple and useful to catch the behav- 
ior of bicxciton lowest level in CuCl even in the nano-to- 
bulk crossover regime, because the exciton and bicxciton 
states in CuCl has been well analyzed by the bipolari- 
ton theory 5 ^ 6 , and RHPS experiments ,5i2L5S While the 
treatment of biexcitons is in general a four-body problem 
with two electrons and two holes and it is usually a hard 
work, owing to the above mentioned simple treatment, 
we can easily discuss the polarization correlation of pho- 
ton pairs emitted from the bicxciton lowest level, which 
has no angular momentum. 

Moreover, by the use of the dyadic Green's function 
for the wave equation of electric field, we can consider 
the surroundings of excitonic material, such as an opti- 
cal cavity consisting of two DBRs. In order to extract 
the scattering fields, instead of using the input-output 
relation , 41 ' 59 " — we consider the definition of Green's 
function and commutation relations of fluctuation opera- 
tors. This simple treatment is valid at least in multilayer 
systems and useful to consider complicated structures. 

In the following subsections, we show our theoretical 
framework to calculate the signal and noise intensities by 
RHPS. We show the Hamiltonian in Sec. Ill Al and the 
equations of motion are derived in Sec. HTbI In order to 
discuss the RHPS, we use some approximations, which 
are explained in Sec. Ill CI The model of biexcitons are 
shown in Sec. Ill Dl In order to solve the equations of mo- 
tion, we use the Green's function technique explained in 
Sec. Ill El Finally, we derive the expression of observablcs 
in Sec. iHFl 



A. Hamiltonian 

Our theoretical framework is based on the QED theory 
of excitons. 44 The Hamiltonian is written as 



H re 



(1) 



where H ex describes the excitonic system, H rcs repre- 
sents a reservoir for the nonradiative damping of exci- 
tons, iJj n t is the exciton-photon interaction, and H em 
describes the electromagnetic fields and a background 
dielectric medium as discussed in Ref. H3 and also used 
in Ref. H3. In order to discuss the biexciton-associated 
RHPS, we consider an exciton-exciton interaction with 
coefficient V^u-y y ■ Namely, the Hamiltonian of exci- 
tonic system is written as 

H eK = Y,^J>lK + \ Yj V^, v <blbtb v \,, (2) 



where is the annihilation operator of an exciton in 
state fj, and is its eigenfrequency. We treat the exci- 
tons as pure bosons satisfying 



0, 



(3a) 
(3b) 



and their non-bosonic behavior is described by the 
exciton-exciton interaction, the second term in Eq. 
The reservoir H rcs is written as 



E 



g ^n)d^n) + g ;(Q)dl(o)\}, (4) 



where d^(f2) is the annihilation operator of harmonic os- 
cillator with frequency fl interacting with excitons in 
state n, and g^(f^) is the coupling coefficient. The os- 
cillators are independent with each other and satisfy the 
following commutation relations: 

[<f M (fl), 4 (n')} = 6^,6(f2 - tf), (5a) 
[d ll (O),d fl ,(O')} = 0. (5b) 

Further, H m t is simply written as a product of electric 
field E(r) and excitonic polarization P cx (r): 



H int = - I dr P ex (r)-E(r). (6) 
Here, the excitonic polarization is represented as 



(7) 



where the coefficient V ^(r) is expressed by the exciton 
center-of-mass wavefunction g c ^{r) and unit vector e M of 
polarization direction as 

7> M (r) = Me Mfl «(r). (8) 

The absolute value of M can be evaluated by the 
longitudinal-transverse (FT) splitting energy Z^lt = 
|M| 2 /e ^bg of excitons, the vacuum permittivity e - an d 
the background dielectric constant e^g °f the excitonic 
medium. 



B. Equations of motion 

According to Ref. H or the QED theories of dispersive 
and absorbing media , 40 ' 41 ' 64 the equation of motion of 
electric field E(r) is derived in frequency domain as 



V x V x E + (r,u) 



-e hg (r,uj)E + (r,uj) 



(9) 
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Here, /zo is the vacuum permeability, £b g (r,cj) is the di- 
electric function of the background medium with arbi- 
trary three dimensional structure. We write an operator 
with a check (") in frequency domain. Jq(t,uj) describes 
the fluctuation of electromagnetic helds and satisfies 



[Jo(r,w),{J (r>'*)}t] = [Jo(r, W ), J (r',-n/)] 



8{ui - u')5(r - r') 



Im[e bg (r,w)] 1 



(10) 



In the same manner as in Rcf. |4J, we obtain the equation 
of excitons' motion in frequency domain as 

[Hujfj, — tvjj- i7cx/2] bn(u) 

dr V;(r)-E+(r,u 1 )+V l ,(u J ) 

/oo iut 
dt bt(t)b v ,(t)K'(t), 
-oo ^ 



fl ' ,v 



(11) 



where 7 0X is the nonradiative damping width defined in 
terms of {g^fi)} as shown in Eq. (D7) of Ref. HH, and 
T> l j l {ui) represents the fluctuation by the damping satisfy- 
ing 



[iy W ),{2V(a/*)}t] = [©»,£„,(-£«/)] 



= <S M , M '<5(w - u> ) 



2vr 



(12) 



The last term on the right hand side of Eq. ((TTj) is the 
nonlinear term due to the exciton-exciton interaction. 
Here, we define a new operator 



(13) 



which annihilates a biexciton (excitonic molecule) in 
state A and describes a two-exciton eigen state i?jjg) by 
applying it to the ground state |g) of matter system. The 
coefficient F\ tlltV is invariant by the exchange of two ex- 
citon indices as 



Further, it is ortho-normal 



(14) 



(15) 



/ ( . V 



and also has a completeness 



(16) 



From the excitonic Hamiltonian _ff ox [Eq. @], the coef- 
ficient F\ ^ M and cigen frequency f2\ of biexciton eigen 
state A should satisfy 

(hoj^ + hw^Fx.^.y + ^2 V^-^yFx^'y = hQ\F\^ v . 

H',u' 

(17) 

By using Eqs. (fT4"]l and ([T^|. we can rewrite Eq. (|13|) as 



(18) 



Therefore, from this relation and Eq. (|17[) . we can rewrite 
Eq. JTl]) as 



[7k^ - fvjj - i7ex/2] 6 M (w) = / drV*Jr)-E + (r,u)+I> li ( : 



/OO 
do/ {S^w'-w)}^"') 



(19) 



On the other hand, by deriving the equation of motion 
for B\ and by using the above relations, we get 

/OO 
da/ (hw v - Hu)')b v (u)')bp(w - u'). 
-oo 

(20) 

In principle, the biexciton RHPS process is described by 
the three equations of motion ©, and (f2"0)l , and 

commutation relations (fTU|) and (|T2"j) . However, in the 
actual calculation, we use the following approximation. 



C. Approximation for RHPS process 



We suppose that a coherent light beam resonantly ex- 
cites biexcitons and their amplitude is large enough com- 
pared to the vacuum fluctuation. In this situation, if we 
do not consider the other higher processes, the biexciton 
operator in the nonlinear term of Eq. (|19[) can be replaced 
by its amplitude B\(lj') = (Bx(uj')). Further, we replace 
6„(o/ — lu) in the nonlinear term by bu (w' — to), which 
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satisfies the linear equation 

[hu^ -fvjj- i7c X /2] b^\u>) 

= J drV;(r)-E + (r,u;)+V,(u J ). (21) 

Simultaneously solving this equation and Eq. b^\u>) 
can be expressed by the fluctuation operators Jq(t,<J) 



By solving this equation and Eq. ([§]), we can repre- 
sent E + (r,ui) by the fluctuation operators Jo(r,u>) and 
T> fl (i J j). This calculation is also straightforward by using 
the Green's function. 

For the calculation of B\(to), we suppose that the biex- 
citon amplitude is not decreased by the scattering, be- 
cause its contribution is small compared to the pumping 
light. Under this approximation, by phenomenologically 
introducing a damping constant 7b x , the biexciton am- 
plitude is obtained from Eq. (|20[) as 

/OO 
dw' (hu, - Huj')(bW(uj'))(bW(u> - u% (23) 
- OO 

where can be calculated from Eqs. (j9]) and (|21| 

by considering an incident light beam as a homogeneous 
solution of Eq. ^ . Under the weak bipolariton regime^ 
where the coupling between exciton-polariton and biex- 
citon is small enough compared to their broadening, the 
approximated expression (|23p of biexciton amplitude is 
sufficient for the discussion of RHPS process. While 
Savasta et al. considered the equation of motion of pro- 
jection operators ; 12 ' 39 they also used a similar approxi- 
mation for the treatment of biexcitons under the detailed 
verification of its validity. 



and r D l j L {uS). Its calculation is straightforward by using 
the Green's function technique as will be shown in sec- 
tion III El Under the above approximation, Eq. (fl~9|) is 
rewritten as 



(22) 



It is well known that the lowest level of biexcitons in 
CuCl is singlet and has zero angular momentum, because 
of the exchange interactions between two electrons and 
between two holesAi Since we suppose the resonant two- 
photon excitation of the lowest level, we only consider 
the lowest relative motion of biexciton in our calculation. 
Further, according to the RHPS experiments in Ref. l54l . 
the lowest biexciton state mainly consists of Is excitons, 
and the contribution from the higher exciton states was 
estimated in the order of 10 -4 . Therefore, we consider 
only Is relative motion of excitons, which has a degree of 
freedom of polarization direction £ M = {x, y, z}. Accord- 
ing to the exciton and biexciton states in bulk CuClfii 
the lowest biexciton level \J = 0, M = 0)b x with zero 
angular momentum is represented as 

| J = 0, M = 0) bx = -{|0, 0; 0, 0) 2cx + |1, 1; 1,-1)** 

+ |l,-l;l,l) 2ex - |l,0;l,0) 2o x}, 

(24) 

where \ji, mi; j2, m2)2ex is the two-exciton state repre- 
sented in terms of angular momenta (ji, mi) and m%) 
of two excitons. This expression surely reflects the po- 
larization correlation of photon pairs observed in RHPS 
experiments 1 ^ 3 - and also determines the phase between 
the two states 



[Huj„ - hu> - i7ex/2] S„(w) / dr?i(r)-^(r,w) + ^(w) 



/OO 
da/ {^(u/ -u)}tB A (a/). 
-oo 



D. Model of biexcitons 

Although Fx ^ u and fl\ should be in principle de- 
termined from Eq. (|17p for given nonlinear coefficient 
Vw,n' y , we instead express F\^ u and fi\ by using 
experimental results. This treatment is useful because 
we already know many parameters of the lowest level of 
biexcitons in CuCl by the longstanding experimental and 
theoretical studies^ 



$+ = (\L,R) + \R,L))/V2 (25a) 
= (\H,H) + \V,V))/V2. (25b) 

Here, \L, R) means that one photon is left- and the other 
is right-circularly polarized, and \R,L) is the opposite 
state. \H,H) and \V, V) respectively means that both 
photons are horizontally and vertically polarized. By 
rewriting each exciton state in terms of the polarization 
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direction as 



E. Green's function technique 



|M> = -(|*> + A/2, 

\j,m) = { \l,0) = \z), 

|l,-l) = (|x)-%))/A 



(26) 



Eq. (|24[) is rewritten as 

|J = 0,M = 0) bx 
1 



{|0,0;0,0) 2o x- \x;x} 2e x- \y;y}< 



\z\ z)2ox}, 

(27) 



which also reflects the polarization correlation (|25b[) . 

Considering the relative motion tf'(r) of two excitons 
in the lowest biexciton level, the coefficient is written as 



^,, = Wj dr J dr'V(r')gT(r)g e **(r+r')gr(r), 

(28) 

where g™ (r) and <?; bx (f ) are center-of-mass wavefunctions 
of excitons and biexcitons, respectively, and 



(29) 



represents the polarization selection rule reflecting the 
lowest state of biexciton [Eq. (|27p]. Here, we suppose that 
the Bohr radius of biexciton (1.5 ran in CuCl)^ is much 
smaller than the crystal size, and the relative motion of 
biexcitons is not strongly modified from the bulk one. 
Namely, we approximate the above expression as 



where <P is defined as 



„* J dr g*(r) g e **(r) 



<P = J dr S'(r) 



(30) 



(31) 



|#| 2 represents the effective volume of the lowest biexci- 
ton state. It was estimated by an experiment^ and was 
also used as a parameter in a theoretical worki££ 



Here, we explain how we simultaneously solve the equa- 
tion of motion of electric field [Eq. ©] and that of exci- 
tons [Eq. (JUJ) or Eq. ([22]l]. By using the dyadic Green's 
function satisfying 



VxVx G (r,r',uj) — -£b g (r,u;) G (r,r',u>) = S(i — r') 1 , 



we can rewrite Eq. ((9j) as 



(32) 



E+(r, to) = E+(r, lo)+^u 2 J dr 1 (r, r', u)-P+(r' ,uj), 

(33) 

where E^r^Lo) represents the electric field in the back- 
ground (H em ) system, and it is defined as 



E+(r,uj) = i/xow / dr' V(r,r',w) • Jb(r',w). (34) 



From Eq. ([TD|). E^(r,uj) satisfies^ 



[E+(r,u;),E^(r\u,')] = [E+{r,u>),E+{r',-J)] 



S(cj-uj')^—[G(r,r',Lj)- G *(r,r',w)]. (35) 

l27T 



^ y 

The expression of G(r,r',u>) in planar system (dielec- 
tric multilayer) is already known^i and will be shown in 

Sec. En] 

Substituting Eq. ([33]) into Eq. (|22"]). we obtain the si- 
multaneous equation set for exciton operators under the 
rotating-wavc approximation (RWA) as 



EV(")VM= / dr- V;(r)-E+(r,oj)+V„(uj) 

/CO 
dJ {&< 1 V-">)} t £A(«')> ( 36 ) 
-OO 



A.;/ 

where the coefficient on the left-hand side is defined as 



S^,{u)=[hw, l -hw-i lcx /2}5^ - ^lj 2 J dr J dr' Vl(r)-^(r,r',u;)-V^(r'). (37) 



The last term of Eq. (|37]) represents the self energy due to the retarded interaction through the electromagnetic 
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fields and the longitudinal Coulomb interaction. Further, 
Eq. (|2Tj) for (oj) in the linear regime is also rewritten 
as 



We separate Eq(v,uj) into an input field -E^ IN (r,aj) 

from r to the matter and an output field -E^~ OUT (r, uj) 
from the matter to r as 



J^-WM^M = / dr ?W ' ^o"(r,w) +2J M (w). ^o"(»-,w) = ^ IN (r,w) + £+ OUT (r,w). 



(44) 



(38) 

This simultaneous linear equation set is solved by cal- 
culating the inverse matrix \N (uj) = [ S (uj)] -1 , and the 
commutation relation of bj_P (uj) is derived in Ref. HH as 



6«M,{6«(c'*)}t 
= 5(w - [W^M - , (39a) 



m\uj),m(-uj' 



Further, Eq. (|36|l is rewritten as 



(39b) 



b^uj) ~b^(uj)+ J2 w^ l (uj)(hu fl , + hw u -hQ x ) 



x F x „, 



duj' {b^(J - W )}tB A ( w '), (40) 



and, by substituting this into Eq. (f3"3"| . the electric field 
involved with RHPS is expressed under the RWA and the 
approximations discussed in Sec. Ill Cl as 

E+ (r,uo)~ E+ (r,uj) + Y J £»(r, uj)b^ (uj) + E NL (r , uj), 

f 

(41) 



where 



e^(r,uj) = fi uj 2 J dr'^(r,r',uj)-V^(r'), (42) 



E^ L (r,uj) = ^2 £ li (r-,uj)W^^(uj)(nhJ li , +t\kJ v - %Q\) 

f-i,f-i' ,A,z^ 



x Fx.u.1 



duo' {^V -uj)YB x (uj'). 



(43) 



F. Input and output fields 

Here, we must pay attention to the electric field 
Eq(t, uj) in the background system, which represents not 
only the field from matter to an observing point r but 
also the field from r to the matter. This means that 
the latter contribution must be removed from Eq. (|4"Tj) 
to calculate observables, while the other terms involving 
(uj) and E^(r,uj) represents the components emit- 
ted from the matter. While such a calculation has been 
usually treated by the input-output relations ; 41 ' 59 " — we 
use the following treatment based on the dyadic Green's 
function G (r,r',ui) for Eq (r,ui). 



By considering the causality, the output field at time t 
should have a correlation only with fields at t' < t, and 
the commutation relation should be written as 



poo poo 



icut — \oj' t' 



E Q (r,t),E (r',t' 



8(uj — uj )- 



i2n 



G (r,r',uj), 



(45) 



where we use the fact that the dyadic Green's function 
G (r, r', uj) satisfying Eq. (|32|) is the retarded correlation 
function of the electric field: 44,68 



- /i ^ 2 G" (r,r',uj) 



(46) 



In the same manner, the input field at t should have a 
correlation only with fields at t' > t, and the commuta- 
tion relation is derived as 



E+ m (r,uj),Eo(r',u>') 

^ poo ft' 



(2tt)' 



dt' J 



iujt— iuj't' 



E (r,t),E (r',t') 



= -5(uj-uj')^—C*(r,r',uj 



i2tt 



(47) 



Actually, Eqs. ([45J and fl47|) reproduces Eq. (|35|). By 
using the output field Eq OVT (t, uj) in the background 
system, we define the scattering field excluding the input 
one as 

^HHPs( r ' UJ ) = E + (r,uj) - E+ m (r,uj) 

= E+ m (r,uj) + E+ L (r,uj), (48) 

where E^ m (r,uj) is the linear component of the electric 



field excluding the input field as 

,(r, uj) = E+ olJT (r, uj) + J2 ^(r, ^ (w). (49) 



E + , 



By deriving commutation relations of £?lin (^j an d 
E NL (r,uj) from Eqs. (|3"5]h (|3"9")l . and (g5]), we can eval- 
uate the observables of RHPS. 



III. PRACTICAL CALCULATION 

Next, we apply the theoretical framework discussed in 
the previous section into multilayer systems embedding 
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FIG. 2: (a) Center-of-mass wavefunctions of excitons and 
biexcitons in a CuCI film. Simple sinusoidal functions vanish- 
ing at surfaces are supposed, (b) Cavity structure considered 
in Figs. [9] and [TO] The Bragg mirrors consist of PbBr2 and 
PbF2. On the transmission side, a high reflectance is achieved 
by a mirror with 16 periods to suppress the leakage of pho- 
tons in this direction. On the incident side, only 4 periods 
are supposed to guarantee rapid radiative decay of entangled 
photons. 



a CuCI layer, and derive expressions of one-photon scat- 
tering intensity and two-photon coincidence intensity by 
RHPS. An incident light beam propagates along z axis 
(perpendicular to the surface), and photon pairs emitted 
into x — z plane are considered (in-plane vector is in x 
direction). We suppose that center-of-mass motions of 
excitons and biexcitons are confined in the CuCI layer 
with thickness d existing at < z < d. Since we consider 
a large enough thickness d compared to the Bohr radii 
of exciton (0.7 nm) and biexciton (1.5 nm)^ the rela- 
tive motions of excitons and biexcitons are not strongly 
modified from the ones in bulk crystals, and all the infor- 
mation of the relative motions are described by factors 
M and <P in Eqs. © and flU). As seen in Fig. EJa), the 
center-of-mass wavefunctions of excitons and biexcitons 
are expanded by a series of sinusoidal functions as 



(r) 



CX 

9k.rr 



(r) 



e*™* 2 
0(*)-^Y^sin(g m 2), 



where 9(z) is unity for < z < d and zero otherwise, 
k is the in-plane wavenumber, S is the normalization 
area in x — y plane, and q m = rrnr/d is the confinement 
wavenumber in z direction for m = 1, 2, .... We con- 
sider eh g (r,uj) as a discontinuous step-like function in 
z direction representing the background dielectric con- 
stant in each layer, and it does not depend on u nor t*|| . 
In the case of multilayer structure in Fig.[^b), Eb g (r,ui) 
gives the background dielectric constant e£g for excitons 
at CuCI layer, otherwise it gives the dielectric constant 
of each layer. According to Ref. |63, if z' is in the CuCI 
layer, the dyadic Green's function satisfying Eq. ([3"2"| is 



expressed as 



G - k (z,z ,w) 



g lk(_X X ) ^ y 

drn / drd G(r,r',oj) 



II 



S 



'i9k cx 
lZk bg 



-^Sh<^-A (51) 

where fc£* 2 = s^ui 2 /c 2 — k 2 and is the unit vector in 
£ direction. When z is in layer j with dielectric constant 
Si, the tensors in Eq. (fBTj) are written as 



g\(z,z',u) = e y e y Q^(z,z',u), 
%%(z,z',cj) = 



D Z D' Z iD z V 







ifcD' o P / VW^ 2 l c2 



(52) 



(53) 



where B z = d/dz and B' z = d/dz'. Eqs. {52) and ([53]) re- 
spectively describe the propagation of V- and H-polarizcd 
fields, and, according to Ref. IoTI GX \z, z') is expressed 
as follows. When < z < d is in the CuCI layer, 



,V/H 



(z,z') 



e ik ^R V/H 



L 

e'" b e 



MV/H 

M V ' H . 
(54a) 



When z is in the leftmost semi-infinite region 

r V/H, is 



-ikLZrpV/H 

1 L 



When z is in the rightmost semi-infinite region, 



M V ' H . 
(54b) 



g][ /H (z,z , ) = fP**f i 



R 



M V ' H . 



(54c) 

Here, R^^ represents the generalized reflection coefH- 
(50) cient for V7 ii-polarized field from the CuCI layer against 
the left(right)-hand neighboring, and ^^R) * s ^ ne § en ~ 
eralized transmission coefficient from the CuCI layer 
to the left (right) most region. The derivation of these 
coefficients is shown in Ref. [|37l Further, k^^ 2 = 
£_l(_r)W 2 /c 2 — k 2 is the wavenumber in the left(right)most 
region with dielectric constant £_l(_r), and the factor is 

defined as M V ' H = [1 - R V L /H R V J H ^ d )-\ 

From Eqs. (f50|) and ([5"T]) . we can evaluate the coefficient 
matrix S (ui) [Eq. (|37p] and numerically calculate the 

inverse matrix W(w) = [ S (w)] _1 . From Eq. ([55]). the 
amplitude of excitons is obtained in linear regime by 

v?V)> =£WWM / dr?;,(r).(£+(r,,)). (55) 
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Here, (Eq (r, uj)) represents the amplitude of electric field 
in the background dielectric system H om , and can be de- 
rived by the standard transfer matrix method^ in the 
case of dielectric multilayers. For simplicity, we consider 
a monochromatic laser light with frequency uj- ln with in- 
plane wavcnumber fc; n . Concerning the pump power L m 
((Eq) oc VTin): there is a scaling law for the intensity of 
entangled photons as explained below. In the present pa- 
per, since we only consider the Is exciton and the lowest 
biexciton level, the exciton states arc labeled by polariza- 
tion direction £ = {x,y, z}, in-plane wavcnumber k, and 



index of center-of-mass motion mas/i = {£, k,m}, and 
the biexciton states are labeled by A = {k, m}. Consider- 
ing the conservation of energy and in-plane wavevector, 
the amplitude of biexcitons is evaluated by Eq. (J23J) and 
we write it as 

%m( W ) = k,2hJ(u - 2w m)#2fc ln ,m( 2w in)- (56) 

Further, the linear and nonlinear components of the scat- 
tering field [Eqs. (|4"9")l and (|4"3"1) ] arc simply rewritten as 



(57) 



where the coefficients are evaluated by the following quantities and functions 



£ i,k,mi z ^) = Hou 2 M^ J dz' ^k(z,z',uj) ■ e^sm(q m z')9(z'), 



(58) 



(59) 



^^L.,fc,m( z ' u 'to,w) — ^2 £(" ) k,m"( Z i UJ )^ r {(",k,m"},{e,k,m'}( UJ ) 

X k,m' + ^e,25 ilx -fe,m — ^2fe in ,n)^{2fe in ,n},{?',fc,m'},{C,2fc in -fc,m}^2fc in ,n(2^in), 

F {2k ln ,n},{e,k, m >},{t;,k<, m } = S^'k'^hn-k® (jj) J dz 9 ( z ) sin(g„2:) sin(q m z) wa.(q m ,z), 



(60) 
(61) 



- 5- k - k ,ii u]' 2 \M\ 2 J dz J dz 1 e{z)6{z')e i ■ tf s («, z', u) ■ e e sm(q m z) sm(q m ,z'). (62) 
Further, from the commutation relations and the following relations are derived for uj > and ui' > as 

; (z,z',o;i n ,a;), (63a) 



LN 



(63b) 



S- LS (r, w ), £+ L - k , (r>, = (w - f^O, ^ «), 



where the tensors arc defined as 



V, 



(64a) 
(64b) 



£,£',771,77l' 
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) = ^ Z ^iCfc.mC^) [ W/ {{,2fc in -S I m},{f',2fc ln -S,m'}(2Win - w) 
" M/ {V,2fc i „-fe,m'} : {C,2fc in -fc,m}( 2a; in ~ W ) 



.fc, m '( Z '' W in,w). 



(66) 



From these commutation relations, we calculate the one-photon scattering intensity and the two-photon coincidence 
intensity. When the background field E j, (z) is in vacuum state in the scattering direction determined by k and only 
has the quantum fluctuation, we obtain the following relations for the initial condition |0) without considering the 
perturbation by the exciton-exciton scattering: 



E+(z,u)\0) = E+ t (z lW )|0) = E+(z,u)\0) = E^Jz,w)\0) = 0. 



(67) 



When we measure the one-photon scattering intensity in the direction k at position z with polarization direction £ 
and frequency w by resolution Aw, the intensity is written as 



^ m! ,:.(^Win,u;) 



ui+Aui/2 



ui — Au/2 



Aw 



(68) 



where E T 



>g g , is the £ component of -E^ HPS ^ and [• ■ ■ extracts (£, £') component of the tensor. Here, it is worth 

noting that this one-photon scattering intensity is proportional to ii n 2 , the square of the pump power, reflecting the 
power dependence of the biexciton creation. Further, the z dependence of this function only represents the scattering 
direction to left or right hand side, if the leftmost and rightmost regions are non-absorptive. When we measure the 
two-photon coincidence between the scattering fields of (£i, ki, Zi, wi) and (£2, k 2 , z 2, ^2), the intensity is calculated 
by 



C 



(2) 



uix+Aui/2 



oji-Alj/2 



ui 2 +Aui/2 



u 2 -Auj/2 



dw 2 dw 2 ' 



^ 1 ,?2,fe 1 „,fc 1 ,fc 2 (Zl ' Z2 ' Win ' Wl ' W2) = 

X (-^RHPS.fei,?! W 'l)^RHPS,fc 2 ,« 2 ( Z2 ' W 2)^RHPS,fc 2 ,6 ^ 2 ' W 2 )^RHPS,fci ,6 fa'^l))- 

By using the above commutation relations, we finally get 



(69) 



C 



(2) 



l^i.Jfel^ 1 '* 2 '^'" 1 '" 2 ) 

+ (Au>)% ir J(oj 1 ,u 2 )\^- ki (z 1 ,Z2,u 1 )] \tf^- k2 (z 2 , Zl ,u> 2 ) . 

L ' J «1>C2 L ' J 6,S1 



(70) 



Here, the function 8(lo, w') gives unity for w ~ w' and zero 
otherwise. The first term represents the signal intensity, 
i.e., the number of correlated photon pairs, which satisfies 
the energy conservation u)\ +uj 2 — 2wj n by resolution Aw, 
and the intensity is calculated as 



/-i(2)S / \ 



(Aw) 2 



LN 



fe„,2fc in -fei( z 2,-Zi,Win,2a;i n - wi) 



6,6 



(71) 



This expression is invariant for swapping the two observ- 
ing conditions, and it is also proportional to Ii n 2 , because 
an entangled-photon pair is emitted from a biexciton. On 
the other hand, the second term in Eq. (|T0"|) has a finite 
value for arbitrary pair of wi and oj 2 , and represents the 



accidental coincidence of emitted photons from indepen- 
dent two biexcitons, because this is just the product of 
two one-photon scattering intensities as 



6 



uh n M (Zl ' Win ' ^ C ilh n M ( Z 2^in, W2) (72) 

This is also invariant for swapping the two observing 
conditions, and proportional to /; n 4 . The third term in 
Eq. ([70]) represents the interference between the two ob- 
serving point, and has a value only for wi ~ uj 2 . There- 
fore, we neglect this term in the following discussion. 

According to Sec. 3.10 in Ref . HT1. we suppose the trans- 
lational masses of excitons and biexcitons are, respec- 
tively, m cx = 2.3too and Wbx = 2.3m cx , where mo is 
the free electron mass. These masses were measured by 
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RHPS experiments! 11 ' 50 ^ 1 However, in our calculation, 
we do not consider the difference of the mass of longitu- 
dinal excitons from that of transverse one. From Sec. 3.2 
in Ref. [Til the transverse exciton energy at band edge is 
fiwT = 3.2022 eV, LT splitting energy is Z\lt = 5.7 meV, 
and background dielectric constant of CuCl is e^g = 5.59. 
Further, according to Sec. 3.7 in Ref. [TTl the binding en- 
ergy of biexciton lowest level is Zibx = 32.2 meV. The 
energy of excitons including the center-of-mass kinetic 
energy is written as 

h 2 - 

^,k,m = + o ( fc2 + I™) ■ (73) 

The energy of biexciton is 

h 2 

m %m = 2hui T - ^bx + (k 2 + q m 2 ) . (74) 

We use the other biexciton parameters reported in 
Ref. [53l : The phenomenological damping width is 7b x = 
h/50 ps = 13.2 fieV, and the effective volume is |^| 2 = 
(4000/2) x (0.541 nm) 3 /4 = 80 nm 3 , where 0.541 nm is 
the lattice constant of CuCl, and 4000 is a parameter 
representing the nonlinear strength. In most of all cal- 
culations, we consider the exciton nonradiative damping 
width as 7 CX = 0.5 meV. 

Because of the translational symmetry in x — y plane, 
the in-plane wavenumber of the system is conserved. In 
the following discussion, we suppose that the pump field 
is perpendicular to the layers, and biexcitons have zero 
in-plane wavenumber. Then, a scattered photon with k 
makes a pair with the one having — k. However, their 
frequencies arc different in general satisfying the energy 
conservation lo\ + LJ2 = 2u>- m . In the present paper, we 
define the scattering angle 9 as fc = (ojt/c) sin 9, which is 
approximately equal to the scattering angle in vacuum. 

IV. RESULTS 

By using the theoretical framework discussed in the 
previous sections, we calculate the scattering spectra by 
bulk crystal and by thin film in Sec. ITVAl In Sec. HVTfl 
we discuss the difference of entangled photon generation 
by thin film from that by bulk crystal, and show the 
thickness dependence of generation efficiency and perfor- 
mance by RHPS in CuCl. Finally, we discuss the the 
generation from a DBR cavity embedding a CuCl layer 
in Sec. ITVCl 



A. Scattering spectra 

Fig. [3] shows forward (transmission side) scattering 
spectra of RHPS from a CuCl film with thickness d = 
7 am. We plot C, 1 } „ E , , % . a (z > d,u>i n ,uj) as a 
function of tu for scattering angles 8 = 0°, 30°, and 60°, 
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FIG. 3: Forward scattering spectra from a CuCl film with 
thickness d = 7 fim. The film exists in vacuum, and the 
pump beam is perpendicular to it. The pump frequency Wi n 
corresponds to the two-photon absorption due to the biexci- 
ton. The results for scattering angles 9 = 0°, 30°, and 60° 
are plotted with different lines as functions of scattering fre- 
quency lj. 



and the spectra are summed for all the polarization di- 
rection £ = {x,y,z}. The CuCl film exists in vacuum, 
and the pump frequency is tuned to the two-photon ab- 
sorption involving biexcitons as hw m ~ huj^ — Ac- 
tually, hu!- m is not exactly Swt — A\, x /2, because we must 
also consider the phase-matching condition (wavevector 
conservation) between two polaritons and a biexcitonAi 
Since the shapes of scattering spectra do not depend on 
the input power 2j n , we plot the scattering intensity with 
arbitrary units. The decay paths of biexcitons are de- 
picted in Fig. [J] As seen in Fig. [31 at 9 = 0°, we have 
multiple peaks at fku — huj^ — — zibx/2 = —16.1 meV 
and a single peak at fkj — Hu>t — — <4bx = 32.2 meV. The 
latter is called Mt peak, which is emitted by the biex- 
citon relaxation into transverse exciton level (exciton- 
like polariton) The remaining polariton with frequency 
lu ~ wt propagates backward, but it cannot go outside 
the film because of the absorption. On the other hand, 
the multiple peaks at Hui — Swt — —16.1 meV originate 
from the biexciton relaxation into two polaritons, and 
the peak structure is due to the interference inside the 
film with d = 7 /xm. Increasing the scattering angle, the 
entangled peaks are split into lower and higher energy 
sides satisfying the energy and wavevector conservations 
as discussed in Ref. These peaks are the LEP and 
HEP, and the intensity of HEP is usually smaller than 
that of LEP, because of the strong absorption near the 
bare exciton energy wt- The angle dependence of the 
peak positions obeys the discussion of Ref. [37]. The peak 
at hoj — Hojt — — ^bx — ^lt is called Ml, which is emit- 
ted by the biexciton relaxation into longitudinal exciton 
state. The emitted photon cannot go outside when 9 = 0° 
because it is polarized in z direction (longitudinal), and 
the remaining exciton also cannot go outside due to the 
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FIG. 4: The scattering spectra of horizontal(H)- and 
vertical(V)-polarizations are shown. The film thickness is 
d = 7 /mi, and the scattering angle is 6 = 60°. The other 
parameters are the same as in Fig. [3] 



strong absorption even for 6 > 0°. 

Fig. |4] shows the polarization-resolved scattering spec- 
tra. The film thickness is also d = 7 /im and the scat- 
tering angle is 9 = 60°. H and V represent the horizon- 
tal and vertical polarizations, respectively, with respect 
to the scattering plane. The Ml peak consists of only 
H-polarized light, because V-polarization does not con- 
tain the longitudinal component. Concerning the LEP 
and HEP peaks, these scattering intensities depend on 
the polarization. We generally get this behavior at a 
non-zero scattering angle, because the reflectance at the 
surface is in general different for the two polarizations. 
When we resolve the spectra with circular polarizations, 
the spectra of left- and right-polarizations are the same 
for any scattering angles and frequencies. 

Fig. [UJa) shows the polarization-resolved scattering 
spectra for thickness d = 200 nm. The scattering an- 
gle is 6 = 60° and the pump frequency is tuned to excite 
the m = 6 biexciton state. Compared to the spectra for 
bulk crystal in Fig. £fj there are more than two peaks 
in the LEP-HEP frequency region. The peak positions 
are different for H- and V-polarizations, and they do not 
obey the angle-frequency relation for bulk crystal^ The 
spectral shape can be interpreted by the exciton-photon 
coupled mode a 17 ' 18 ' 21 ' 32 ' 33 ' 35 in the thin film, which have 
been discussed in relation with the radiative decay of ex- 
citons in nano-to-bulk crossover regime^ - — Due to the 
breaking of translational symmetry in z direction, the 
lower and upper polaritons in bulk material are no longer 
good eigen states, but instead we obtain the exciton- 
photon coupled modes with discrete energy levels and 
radiative decay rates in the case of thin films. A cre- 
ated biexciton spontaneously decays into these coupled 
modes with emitting a photon conserving the energy and 
in-plane wavevector. By using the method in Ref. [l~8l . 
we calculated the exciton-photon coupled modes with V- 
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FIG. 5: (a) Polarization-resolved scattering spectra for thick- 
ness d = 200 nm and scattering angle — 60°. The film 
exists in vacuum and the pump frequency is tuned to reso- 
nantly excite the m — 6 biexciton state, (b) Dashed lines 
are the dispersion relation of exciton-polariton in bulk crys- 
tal. The horizontal bars represent the exciton-photon cou- 
pled modes with V-polarization in the film with d — 200 nm. 
The bar length is the sum of radiative and nonradiative de- 
cay rates and the center is the resonance frequency, (c) The 
H-polarized modes are plotted. Because of the breaking of 
translational invariance in the z direction and the non-zero 
scattering angle, the longitudinal excitons are also optically 
allowed. 



polarization in the film with d = 200 nm as shown in 
Fig. EJb) and the modes with H-polarization are shown 
in Fig. 03c). The dashed lines represent the dispersion 
relation of exciton-polariton in bulk crystal, and the hor- 
izontal bars are the coupled modes in the thin film. The 
length of each bar represents the sum of radiative and 
nonradiative decay rates, and the center is the resonant 
frequency. Since the H-polarized modes includes the 
longitudinal components, there are also the exciton-like 
modes with longitudinal exciton energy. The higher fre- 
quency parts of the scattering spectra in Fig.JSJa) appar- 
ently reflect the structure of the coupled modes shown in 
Figs. 03b) and (c), and the peaks in lower frequency part 
appear with satisfying the energy conservation. In this 
way, the scattering spectra of thin films are completely 
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FIG. 6: The two-photon coincidence (signal) intensity is plot- 
ted as a function of scattering frequency cj of one photon. 
The spectra are resolved by the polarization of the two pho- 
tons. CuCl films with thicknesses d = 7 /im and 200 nm are 
considered, and the scattering angle is — 60°. The same 
pump power is supposed for both the thicknesses, (a) The 
polarization directions of two photons are resolved in H- and 
V-axes, and HH and VV pairs are correlated. The pairs with 
HV and VH polarizations are not generated from the lowest 
biexciton level in a film, (b) The polarization directions are 
resolved for circular polarization base. Not only the pairs of 
left (L) and right (R) circularly polarizations are obtained, 
but LL and RR pairs are also generated for 9 > 0. The spec- 
tra of LR and RL are the same, and those of RR and LL are 
also the same. The spectra of 200 nm are magnified by factor 
80 in both (a) and (b). Since the spectra are symmetric about 
ojin, only the higher frequency side uj > cj- m is plotted. 



different from the bulk one, and they depend on the film 
thickness, surroundings, and in-plane wavenumber obey- 
ing the change of exciton-photon coupled modes as dis- 
cussed in Ref. [l8l . Furthermore, in contrast to the spectra 
for bulk crystal in Fig. [4) the emission near the exciton 
resonance uj — wt can go outside the film, because of the 
large radiative decay rate in the thin film. From these 
results, the measurement of scattering spectra of RHPS 
can be considered as a powerful tool 3 -^ to observe the 
exciton-photon coupled modes in nano-structured mate- 
rials in addition to the previously performed nonlinear 
optical responses ! 17 ' 33 



B. Entangled-photon generation 



Next, we discuss the entangled-photon generation by 
RHPS process. Fig. [S] shows polarization-resolved spec- 
tra of two-photon coincidence measurement. We plot 



only the signal intensity C 



(2)S _ 

£i 7^2, &in— 3 fei = (wt/c) sin 9 



(zi > 

d,Z2 > d, cji n ,w) as a function of scattering frequency to 
of one photon (the other photon has a frequency 2io- m — lo). 
In Fig.^a), the pairs with £i = £2 = H and £1 = £2 = 
are considered, and the pairs with different polarizations, 
such as HV and VH, have no correlation, because the 
lowest biexciton level with zero angular momentum is ex- 
cited. Two film thicknesses 7 /mi and 200 nm are consid- 
ered, and the parameters are the same as in Figs. 2] and 
respectively. In the two calculation, we considered the 
same pump power and the ratio of signal intensities 
do not depend on /; n . Since the spectra are symmetrical 
about Wi n , we show only the higher energy part co > Lu- nl . 
As similar as the scattering spectra in Fig. [4j the intensi- 
ties of HH and of VV are not the same in general in the 
case of non-zero scattering angle. Therefore, the ideal 
entanglement in Eq. (|25b[) is not generally obtained, and 
the entangled state also have RR and LL components, 
whose spectra are shown in Fig. [HJb) . Although this is 
a general property of bulk crystals, the situation is dif- 
ferent in the case of nano-to-bulk crossover regime. As 
seen in Fig. Eta), we can obtain the same signal intensi- 
ties for HH and VV at frequencies hu — huT = —9.8 and 
—7.7 meV for d = 200 nm, and the signal intensities of 
RR and LL become nearly zero at frequency —8.8 meV 
in Fig.[6tb), while it is slightly different from the peak fre- 
quency —8.2 meV of LR spectrum. These results mean 
that the state of emitted photon pairs can be modified 
by tuning the film thickness, scattering angle, and scat- 
tering frequency in the case of thin films. For example, 
at frequency —8.8 meV for d = 200 nm, we can get the 
entangled state (\LR) + \RL)) / \/2, while the proportions 
of \HH) and \VV) are not equal as seen in Fig. [6ja), 
because the polarization basis of the two photons are dif- 
ferent for 9 7^ 0. On the other hand, at frequencies —9.8 
and —7.7 meV, we get the entangled pairs with the same 
HH and VV proportions, while they contains RR and 
LL components. 

Furthermore, even if the scattering angle is 9 = 0°, in 
contrast to the bulk case, the scattering peaks are not at 
lu = lot in general in the case of thin films. Therefore, 
the maximally entangled photon pairs are obtained by 
the frequency filtering under the observation at 9 = 0°. 
Fig. [7Ja) shows the spectrum of signal intensity (gener- 
ation efficiency) obtained by a CuCl film with thickness 
d = 200 nm for scattering angle = 0°. The proportions 
of HH and VV are the same, and RL and LR paris 
are not emitted. As seen in Fig. E£a), the peaks appear 
not at lo = Lo- m but close to the resonance frequencies 
of the exciton-photon coupled modes shown in Fig. (TJc) 
(but not just at the resonance frequency because we get 
weaker absorption at frequency far from wt). In this way, 
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FIG. 7: (a) The generation efficiency (signal intensity) from 
a CuCl film with thickness d = 200 nm is plotted as a func- 
tion of scattering frequency w. The biexciton state m = 6 is 
resonantly excited, and the scattering angle is 9 = 0°. (b) 
The spectrum of the corresponding performance P. The gen- 
eration efficiency and the performance are normalized in the 
same manner as in Fig. [5] (c) The exciton-photon coupled 
modes for 6 = 0° in the film are plotted, while 9 = 60° in 

Fig.ETb). 



compared to the bulk crysta fa 10 ! 13 and also to the sim- 
ple quantum dots^^ the nano-to-bulk crossover regime 
has a variety of degrees of freedom to tune the generated 
state. 

For the high-power generation of the entangled pho- 
tons, the important factors are the generation efficiency 
and also the statistical accuracy, i.e., the amount of un- 
entangled pairs. For a pumping beam with power 7; n , the 
signal intensity S oc C^ s (amount of entangled pairs) is 
proportional to ii n 2 , while the noise intensity N oc C^ N 
(amount of uncntanglcd pairs) is proportional to ii n 4 , 
because an unentangled pair is involved with two biexci- 
tons. This implies that, by increasing the pump power 
Ji n , the S/N ratio decreases in contrast to the increase in 
S'ji 3 - To evaluate the material potential for the generation 
of strong and qualified entangled-photon beams, we intro- 
duce another measure termed "performance" P defined 
as the signal intensity S under a certain S/N ratio a (I{ n 
is tuned to give this ratio). This quantity P = S 2 /aN 
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FIG. 8: Thickness dependences of (a) generation efficiency 
S/Iin 2 and (b) performance P. The scattering angles are 9 = 
0° and the frequencies are W1/2 = win ± + . To suppress the 
interference effect, outside medium is the dielectrics with e^g. 
The results for 7 CX = 0, 0.1, 0.5, and 1.0 meV are plotted 
with different lines. The performance is normalized to the 
ideal quantity. 



does not depend on /;„ and reflects the material potential. 
Fig. [TJb) shows the spectrum of the performance. As 
comparing with Fig.[7][c), the spectrum of P basically re- 
flects the exciton-photon couples modes. However, since 
the spectrum of the noise (twice the scattering intensity) 
is different from the signal one, the Figs. Eta) and Efb) 
are slightly different. The most significant difference is 
the spectra around ui = ui- m ~ ojt — 16.1meV. While 
both S and P mostly reflect the resonance frequency of 
the exciton-photon coupled modes, the performance is 
maximized at uj — uj- m , because P is strongly affected 
by the nonradiative damping, which becomes smaller at 
that frequency. 

Fig. |8] shows the thickness dependences of (a) genera- 
tion efficiency S/L ln 2 and (b) performance P. The shapes 
of P curves do not depend on a chosen a, and the max- 
imum value is normalized to unity. We also plot the 
generation efficiency with arbitrary units, because the 
estimation of absolute signal intensities are sensitive to 
the change of measurement conditions, while the spectral 
shape and thickness dependence do not depend on such 
conditions. For simplicity, we assume that the two scat- 
tering fields are forward and perpendicular to the film 
(6 = 0°) and the frequencies are ui\/2 = w; n ± + . The 
pump frequency is tuned to the two-photon absorption in 
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bulk material. The results for nonradiative decay rates 
7 CX = 0, 0.1, 0.5, and 1.0 meV are plotted with differ- 
ent lines. In order to suppress the oscillation due to the 
interference as seen in Fig. [31 we suppose that the CuCl 
film exists in a dielectric medium with £|*. The results 
for the film in vacuum will be shown in Fig. [9] In Fig. |8l 
the oscillating behavior in the nanometer thickness range 
is due to the biexciton confinement and the modification 
of the energy structure of exciton-photon coupled modes. 
The RHPS effectively occurs when the resonance energy 
of the coupled mode is just equal to half the biexciton en- 
ergy. The maximum performance shown in Fig. |8£b) is 
the ideal quantity, and it only depends on measurement 
conditions, such as resolutions of angle and frequency, 
but not on material parameters. 

As seen in Fig. Eta) and also in Ref. HH, the optimal 
thickness for generation efficiency depends on 7 GX , and 
it is on the order of micrometers or more for CuCl crys- 
tals. However, as seen in Fig. H^b), the performance sig- 
nificantly decreases from the ideal value at a thickness 
of micrometers for nonzero 7 ex , because the nonradia- 
tive decay easily increases the amount of unentangled 
pairs. Therefore, when we use bulk crystals, the gen- 
eration efficiency (generation probability for one pump 
pulse) is limited by a desired statistical accuracy (S/N 
ratio) fi^. However, at a thickness from 50 to 1000 nm, 
as expected, a nearly ideal performance can be obtained 
at particular thicknesses even if 7 ex is nonzero. This is 
because the radiative decay is dominant owing to the 
exciton superradiance) 17 ' 18 and the entangled pairs can 
go outside the film without decreasing the amplitude. 
Therefore, thin films generally show a high performance, 
and this rapid decay is also desired for the high-repetition 
excitation, which also recovers the signal intensity while 
maintaining the S/N ratioJ^ 



C. With DBR cavity 

Although a good performance is obtained at a thick- 
ness of hundreds nanometers, the generation efficiency 
of such thin films is much lower than that of bulk crys- 
tals as seen in Fig. [5{a), and a strong pump power is 
required to achieve a sufficient signal intensity at such 
thickness range. While the superradiant excitons main- 
tain the large nonlinearity (excitonic component)^ this 
low efficiency simply comes from the small thickness (in- 
teraction volume). This problem can be overcome by 
using an optical cavity in the strong-coupling regime, be- 
cause we can control both the interaction volume and ra- 
diative decay rate using two parameters: quality factor 
(Q-factor) of cavity and thickness of CuCl. This aspect 
is different from the simple semiconductor microcavity, 
in which the interaction volume and radiative decay rate 
are respectively enhanced in strong- and weak-coupling 
regimes. 

Although a high generation efficiency can be achieved 
by using a high-Q cavity, we have to simultaneously re- 
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FIG. 9: Thickness dependences of (a) generation efficiency 
S/Iin 2 and (b) performance P. The black lines represent the 
results for a bare CuCl film existing in vacuum. The scatter- 
ing frequencies are the same as that of pumping beam, and 
scattering is forward with 8 — 0°. The gray lines represent 
the results for DBR cavity embedding a CuCl layer shown in 
Fig. H£b). The mode frequency of the optical cavity is tuned 
to ojt, and the scattering is backward with 6 = 180°. In both 
cases, 7ex = 0.5 meV. The performance is normalized in the 
same manner as in Fig. [8] 



alize a rapid radiative decay of entangled photons inside 
the cavity. Therefore, we consider a low-Q cavity as re- 
ported in Ref. [6^, namely, a CuCl layer in a DBR cav- 
ity composed of PbF2 and PhBr2 as seen in Fig. [UJb). 
Here, since the translational invariance is broken at the 
thickness of nanometers, the generated photons can go 
forward and also backward in contrast to the bulk case. 
Therefore, we suppose a high reflectance on the transmis- 
sion side to suppress the leakage of entangled photons, 
and we focus on the backward emission. The DBR cav- 
ity is considered by the background dielectric function 
eb g (r,Lu) in Eq. (|9]). The refractive indexes of PbF 2 and 
PbBr2 are 1.86 and 2.95, respectively. The gray lines 
in Fig. [S] represent the backward emission from the cav- 
ity structure, where the cavity mode frequency is tuned 
to wt, 7ox = 0.5 meV, and the periods of the inci- 
dent and transmitted sides are 4 and 16, respectively 
(Q-factor is 50). This system corresponds to the weak 
bipolariton regim o 15 ' 16 (but the strong-coupling regime 
of excitons and photons), where the energy splitting be- 
tween polariton and biexciton levels is small compared 
with their broadening. This situation is in contrast to 
that in Refs. [ID and \WL where the strong enhancement 
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of cntanglcd-photon generation from a quantum well in 
a high-Q cavity has been discussed on the basis of the 
biexcitonic cavity-QED picture or the strong bipolariton 
picture. As shown in Fig. Eta), since biexcitons are ef- 
fectively created, the generation efficiency is significantly 
enhanced at a thickness of nanometers, and it is larger 
than the maximum value in the case of bare CuCl film 
existing in vacuum (black line). The enhancement also 
occurs when the polariton energy (exciton-photon cou- 
pled mode) is equal to half the biexciton energy, which is 
consistent with the results in Refs. [Tol and [l6l Compared 
with FiglUa), the period of the oscillation is doubled, 
because the RHPS involving biexcitons with odd-parity 
center-of-mass motion is forbidden in the one-sided op- 
tical cavity. On the other hand, as shown in Fig. EJb), 
at a thickness of micrometers, the performance is sup- 
pressed compared with that of the bare CuCl film. This 
is because of the multiple reflections inside the cavity, 
and the scattered fields non-radiatively decay during the 
propagation. In contrast, at the thickness of nanometers, 
particularly at 80 nm, the performance almost maintains 
the ideal quantity. This is due to the enhancement of 
the radiative decay rate by exciton supcrradiance, and 
the enhancement of generation efficiency is simultane- 
ously obtained by the cavity effect in the strong-coupling 
regime. 

Finally, in Figs. [TUta - ) and (b), we show the spectra 
of generation efficiency and performance, respectively, in 
the case of CuCl film with thickness d = 72 nm embed- 
ded in the DBR cavity discussed above. This thickness is 
chosen to achieve a high generation efficiency, while the 
performance at uj = u>- m is 0.86, which is smaller than 
the maximum value in Fig. EJb). However, compared 
to the thin film without the cavity, the generation effi- 
ciency is significantly increased and the high performance 
is successfully maintained due to the rapid radiative de- 
cay. Further, while the pump frequency w; n is assumed 
to the two-photon absorption frequency in bulk CuCl in 
Figs. [TUJa) and (b), we have numerically checked that, 
when uj- ln is correctly tuned to the eigen frequency of 
a confined biexciton mode, the generation efficiency is 
enhanced more dramatically with maintaining the high 
performance in the case with DBR cavity. 

In Fig. [TUtc). the exciton-photon coupled modes are 
plotted with horizontal bars, and we can find that one 
mode with high radiative decay rate exists close to the 
pump frequency w; n (two-photon absorption frequency 
of biexcitons) . This mode corresponds to the lower cav- 
ity polariton, and the strong electric field in the cavity 
enhances the generation efficiency of biexcitons due to 
the cavity-induced double resonance^ Furthermore, the 
generated entangled excitons rapidly decay into photons 
through this polariton mode, which ensures the high per- 
formance at the same time. 

In this way, by using an optical cavity embedding a 
CuCl layer with a thickness of nanometers, we can obtain 
high efficiency and performance simultaneously. To avoid 
the leakage of generated photons, the reflectance on the 




( Pico - h(B T ) / meV 



FIG. 10: (a) The generation efficiency (signal intensity) is 
plotted as a function of scattering frequency uj. A CuCl film 
with thickness d — 72 nm is embedded in the DBR cavity 
considered in Fig. [9j and the other parameters are also the 
same, (b) The spectrum of the corresponding performance P. 
(c) The exciton-photon coupled modes in the film are plotted 
in the same manner as in Fig. Ob). 

transmission side should be high enough, but that on 
the incident side should not be high to obtain a rapid 
radiative decay rate. Once we choose a cavity structure, 
we can numerically determine the optimal thickness of 
CuCl layer, in which an exciton-photon coupled mode has 
half the biexciton frequency, a rapid radiative decay rate, 
and also large exciton component (large nonlinearity) to 
achieve high performance and generation efficiency. Such 
a coupled mode is an unique feature in the nano-to-bulk 
crossover regime. 



V. SUMMARY 

We have developed a theoretical framework for the 
investigation of the biexciton-involved RHPS based on 
the QED theory for excitons 4^ Compared to the previ- 
ous theories) 12 ' 39 our method can be applied to nano-to- 
bulk crossover regime because we explicitly consider the 
center-of-mass motion of exciton. Further, we can dis- 
cuss the polarization correlation of entangled pairs and 
the surroundings of the excitonic layer, such as the DBR 
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cavity structure. While wc considered CuCl films in ac- 
tual calculation, by treating several relative exciton lev- 
els and by correctly calculating the center-of-mass wave- 
functions of excitons confined in finite crystal including 
the effect of dead layer jZl our theoretical framework can 
be applied to other materials in principle. Further, by 
correctly treat the modification of relative motion of ex- 
citons and biexcitons strongly confined in nano crystals 
and also the Pauli's exclusion principle, our framework 
can be extended for the investigation of single quantum 
dot, the deterministic generation of entangled photons. 

We have calculated the scattering spectra of RHPS 
from CuCl films with bulk-like and submicron thick- 
nesses. For the bulk-like thickness, the four peaks called 
Mt, Ml, LEP and HEP are reproduced. On the other 
hand, the scattering spectra of the thin film are sig- 
nificantly modified from the bulk ones, and we found 
that they reflect the exciton-photon coupled modes in 
the thin filmJii&Zi^^ Therefore, the RHPS measure- 
ment is also useful to observe the exciton-photon cou- 
pled modes in nano-structured materials 3 ^ as well as the 
four-wave mixin g 17 ' 29 ' 35 and the two-photon excitation 
scattering^ We also found that semiconductor thin films 
have much degrees of freedom to control generated states 
of entangled photon pairs. 

In addition to the signal intensity of entangled-photon 
generation, we also discuss the performance of the ma- 
terial structure by considering the noise intensity from 
independent two biexcitons. Although the thickness de- 
pendence of signal intensity has a maximum value at par- 
ticular thickness determined by nonradiative decay rate 
of excitons ) 12 ' 39 a high performance is generally obtained 



at thickness of nanometers due to the rapid radiative 
decay of excitons. However, the generation efficiency of 
such thin films is much weaker compared to the bulk one. 
We have demonstrated that, by using a DBR cavity em- 
bedding an excitonic layer in the nano-to-bulk crossover 
regime, the generation efficiency can be enhanced while 
maintaining the high performance. 

For the pursuit of high-power and high-quality but 
probabilistic generation of entangled photons, which is 
essential for the next-generation technologies of fabri- 
cation and chemical reaction^ the biexciton-involved 
RHPS shows a quite high generation efficiency compared 
to that of PDC method. From the viewpoint of the qual- 
ity of generated entangled pairs, the RHPS method can 
show a high performance and a high generation efficiency 
simultaneously by using and an optical cavity embedding 
a CuCl nano-layer. Further, it has much degrees of free- 
dom to control the generated states of entangled photons. 
We believe that our results make a breakthrough in high- 
power and high-quality entangled-photon generation. 
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